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SECTION  1 
INTRODUCTION 

Recently  there  has  been  a  renewed  Interest  In  transmission-line 
theory  and  Its  application  to  the  Internal  Interaction  problems 
Involving  electromagnetic  pulse  (EMP)  excitation  of  aerospace  systems. 

One  new  development  In  this  area  has  been  the  formulation  of  an 
analysis  procedure  to  study  large  interconnected  networks  of  multi¬ 
conductor  transmission  lines.  This  analysis,  which  Is  described  In 
refs.  (1)  and  (2),  and  the  resulting  computer  program  (ref.  3),  will 
permit  not  only  simple  branching  of  transmission  lines  within  the 
network,  but  also  complicated  loop!*-'  of  lines.  Thus,  an  arbitrarily 
Interconnected  set  of  transmission  ,r.es  can  be  analyzed  using  this 
approach. 

The  analysis  of  the  transmission- line  networks  described  In 
refs.  (1)  and  (2)  Is  based  on  the  network  excitation  being  due  to 
lumped  (or  discrete)  voltage  and  current  sources  located  at  a  source 
position  somewhere  along  each  transmission-line  section  (tube).  While 
this  specification  of  sources  may  be  useful  for  certain  applications. 

It  Is  not  particularly  useful  for  EMP  studies,  where  the  transmission- 
line  network  Is  excited  by  an  .ncident,  transient  electromagnetic 
field.  In  the  EMP  case,  not  only  Is  the  transmission- line  excitation 
distributed  along  the  line,  but  the  fundamental  excitation  quantities 
are  the  Incident  electric  and  magnetic  fields  (f  and  5),  not  the 
current  and  voltage  sources.  Thus,  it  Is  necessary  to  modify  the  past 
analysis  to  permit  distributed  field  excitation  of  the  transmission  lines. 


1.  BaiOT,  C.E.,  et  al.,  "N»m<erical  Results  for  Multi conductor  Trans¬ 
mission  Line  Networks,"  AFWl-TR-77-123,  Air  Force  Weapons 
laboratory,  Klrtland  AFB,  NM,  November  1977. 

2.  Baum,  C.E.,  T.K.  Liu  and  F.M.  Tesche,  "On  the  General  Analysis 
of  Multi  conductor  Transmission  Line  Networks,"  AFWL  EMP 
Interaction  Note,  to  be  published. 

3.  Tesr''*  F.M. ,  and  T.K.  Liu,  "User  Manual  and  Code  Description  for 
QV'  A  General  Multiconductor  Transmission  Line  Analysis  Code," 
prr  -ed  for  Air  Force  Weaoons  Laboratory,  Contract  F29G01-78-C-0002, 
l  t  1978. 
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Field  excitation  of  simple  c  en  two-wire  lines  has  been 
considered  by  a  number  of  authors  and  two  sepa»  ite ,  but  equivalent, 
approaches  used.  Taylor,  Satterwhlte  and  Harr  son  (ref.  4)  and 
Smith  (ref.  5)  derive  a  coupling  mode"  based  on  the  Incident  tangen¬ 
tial  electric  fields  on  both  wires  of  the  tr.  emission  line  and  on 
the  short  wires  of  the  loads  at  the  eids  of  the  Ine.  This  approach 
Is  based  on  the  Integral  form  of  Marvell's  equations  as  applied  to 
the  closed  loop  formed  by  the  two  parallel  wires  of  the  transmission 
line  and  the  two  loads  at  the  ends.  In  this  formulation,  there 
appear  distributed  voltage  sources  In  both  wires  of  the  transmission 
line,  as  well  as  voltage  sources  at  both  loads  terminating  the  line. 

A  different  approach  has  been  used  by  Lee  (ref.  6)  to  determine 
the  distributed  field  excitation.  This  Is  based  on  the  differential 
forms  of  Maxwell's  equations  and  yields  distributed  current  and 
voltage  sources  along  the  line,  with  the  voltage  source  being  propor¬ 
tional  to  the  H  field  and  the  current  source  being  related  to  the  f  field. 

Both  of  these  formulations  yield  Identical  results  for  computing 
the  TEM  currents  flowing  on  a  two-wire  line  excited  by  an  Incident 
field.  The  former  approach  has  been  extended  to  the  case  of  multi¬ 
conductor  transmission  lines  by  Paul  (ref.  7)  and  Frankel  (ref.  8), 
and  Is  similar  to  that  discussed  In  this  report.  A  r'lgntly  different 

4.  Taylor,  C.O. ,  R.S.  SatterwMte  and  C.W.  Harrison.  Jr.,  "The  Response 
of  a  Terminated  Two-Wire  Tr  smlsslon  Line  Excited  by  a  Nonuniform 
Electromagnetic  Field,"  -FWL  EMP  Interaction  Notes,  Note  66, 

November  1965;  also,  IEEE,  YransV  A.P. ,  Vol  AP-1$.  pp.  987-989,  1965. 

5.  Smith,  A. A.,  Coupling  of  External  Electromagnetic  Fields  to 
Transmission  lines.  John  Wiley  and  Sons.  New  York,  1$)7. 

6.  Lee,  K.S.H.,  "Balanced  Transmission  Lines  In  External  Fields,"  AFWL 
EMP  Interaction  Notes.  Note  115,  July  1972. 

7.  Paul,  C.R.,  "Frequency  Response  of  Multi  conductor  Transmission  Lines 
Illuminated  by  an  Electric  Field,"  IEEE  Trans.  EMC,  Vol.  EMC-18, 

No.  4,  pp.  183-186,  November  1976. 

8.  Frenkel,  S. ,  Multi conductor  Transmission  Line  Analysis.  Artech  House, 
1977. 
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approach  has  been  employed  by  Kajfez  and  Wilton  In  ref.  (9),  where  the 
concepts  of  reciprocity  have  been  used  to  obtain  the  multi  conductor 
transmission- line  response  to  a  small  aperture  excitation  of  the  line. 
The  method  of  refs.  (4)  and  (5)  has  been  applied  to  multi  conductor 
systems  by  Strawe  (ref.  10),  but  his  report  Is  not  «,.>.ely  distributed. 

The  present  report  discusses  In  detail  the  citation  of 
mul tl conductor  transmission  lines  by  an  Incident  electromagnetic 
field  using  the  differential  formulation.  Section  11  presents  the 
derivation  of  the  equations  describing  the  termite1,  or  load,  current 
responses  of  a  multiconductor  transmission  line.  These  equations  have, 
as  sources,  both  distributed  voltage  and  current  generators  which  are 
induced  by  Incident  magnetic  and  electric  fields.  Section  III  first 
discusses  the  derivation  of  these  local  r.our.es  In  terms  of  the  local 
fields  and  transmission-line  geometry.  The  concept  of  an  "equivalent 
separation"  between  conductors,  as  coBir<nly  used  for  two-wire  lines,  Is 
then  developed  for  an  arbitrary  multiconductor  transmission  line. 
Finally,  In  Section  IV,  the  incident  field  components  which  contribute 
to  the  distributed  sources  are  given  for  an  incident  plane  wave  striking 
the  line  at  an  arbitrary  angle  ->f  incidence. 


9.  Kajfez,  D. ,  and  D.R.  Wilton,  "Small  Aperture  on  a  Multi  conductor 
Transmission  Line  Filled  with  Inhomogeneous  Dielectrics,"  AF0SR- 
76-3025-2,  Air  Force  Office  of  Scientific  Research,  November  1977. 

10.  Strawe,  D.F.,  "Analysis  of  Uniform  Multiwire  Transmission  Lines," 
Boeing  Report  D2-26088-1  under  Contract  F04701-72-C-0210, 

November  1972. 


SECTION  II 

MUITICONOUCTOR  TRANSMISSION-LINE  RESPONSE  TO  DISTRIBUTED  SOURCES 


As  discussed  In  ref.  (I),  the  response  of  a  general  transmission¬ 
line  network  may  be  calculated  by  decomposing  the  currents  on  each 
tube  of  the  transmlslson  line  Into  forward  and  reverse  propagating 
components.  At  every  junction  within  the  network,  a  scattering 
matrix  can  be  derived  to  express  all  scattered  components  of  current 
In  terms  of  the  Incident  components.  These  two  sets  of  relations  can 
be  combined  to  form  a  large  matrix  equation  for  the  Incident  currents. 
This  equation,  called  the  BIT  equation,  can  be  Inverted  numerically 
and  the  Incident  currents  determined.  Through  the  scattering  matrices, 
the  scattered  and,  thus,  the  total  currents  on  the  lines,  can  be 
determined. 

A  basic  element  of  the  above  network  analysis  Is  the  determina¬ 
tion  of  the  propagation  properties  of  the  forward  and  backward  waves 
on  the  line,  as  well  as  their  relative  excitation  by  sources  along  the 
line.  For  the  purpose  of  this  section,  therefore,  we  will  consider 
only  a  single  section  (tube)  of  multi  conductor  transmission  line. 

Consider  a  lossless  section  ot  multi  conductor  transmission  line 

having  no  sources,  as  shown  In  Figure  1.  The  length  of  the  line  Is 

denoted  by  f.  and  It  contains  N  wires  with  the  N+lst  wire  being 

the  reference  conductor.  The  N*J  wires  are  required  to  be  parallel, 

but  not  necessarily  coplanar.  For  ch  a  line,  its  electrical  properties 

are  determined  by  a  caoacltlve  coeff  :ient  matrix,  (C *  )  ,  and  an 

n,m 

Inductive  coefficient  matrix,  ( L ’  _}  ,  which  depend  only  on  line 

n 

geometry  and  dielectric  properties  around  the  line.  For  this  line, 
these  matrices  are  nonsingular  matrices  of  order  N. 

As  discussed  In  ref.  (1),  the  voltages  and  currents  on  this  line 
without  sources  must  obey  a  coupled  set  of  partial  differential 
equations  as 
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( <»„.■>  \  ( 

VC;.«'  \  nn(j's)) 

where  the  notation  (Vn)  represents  an  N-vector  for  the  line  voltage 
and  a  similar  notation  holds  for  the  current.  The  parameter  s  is 
the  complex  frequency  variable,  and  the  tilde  represents  a  Laplace 
transformed  quantity. 

Equation  (1)  can  be  manipulated  into  two  separate  equations 
for  voltage  and  current  vectors.  The  current  equation  becomes 

32(LU.S))  ? 

- ~  5 - s2(c:  )  (L*  )  (I  (z.s))  -  (0n)  (2) 

2  n,m  n,m  n  " 

aZ 


I  <»„(..*» 


which  is  a  one-dimensional  wave  equation  for  the  N-vector  current. 

For  a  lossless  multiconductor  section  imeersed  in  a  uniform, 
homogeneous  dielectric,  the  matrix  product  (C^  „)(L^  m)  1n 
Equation  (21  is  diagonal  and  the  individual  elements  of  the  current 
N-vector  are  themselves  a  solution  to  a  simple  wave  equation: 


(z.s) 


0 


where  v  is  the  velocity  of  wave  propagation  on  the  line. 

A  more  general  line,  however,  does  not  have  a  diagonal  result  for 

the  (C*  1(1*  )  matrix,  although  it  is  possible  to  diagonalize  it 

n  n,m 

through  the  use  of  a  nonsingular  NxN  transformation  matrix,  denoted 
by  (T  )  ,  which  consists  of  the  current  eigenmodes,  (♦  )  ,  as 

n  |W  n  j 

columns.  The  #n's  are  solutions  to  the  eigenvalue  equation 


2 


Un) 

n  j 


(3) 


10 


where 


1 


Is  the  1th  eigenvalue  corresponding  to  the  eigenmode 
By  Introducing  a  change  of  variables  as 


<*«<*•*»  ■  (Tn.«)(V*’s» 


(4) 


where  (1  (z.s))  represents  the  modal  currents,  the  wave  equation 
for  the  modal  currents  becomes 


a2(?n(z.s)) 
~T~ 


dz 


(y_  _)2(1> 
0 1*  0 


(5) 


2  2 
where  (y  .)  Is  a  diagonal  rnatrl'  containing  the  terns  as 
n  |R  i 

elements. 

Since  the  matrix  (•*  ,  :  ^uatlon  (5)  Is  diagonalized, 

•  • 

the  solution  for  the  modal  currents  can  be  expressed  directly  as 
exponential  functions  of  position,  and  the  total  solution  for  the 
line  currents  becomes 


(In(z.s)) 


(T 


n,m 


) 


♦  J* 

(i*)  ♦  e  n*m 


(6) 


where  (v)  and  (a*)  are  N-vectors  which  define  the  amplitudes  of 
n  n 

each  of  the  propagating  modes  on  the  line  and  which  depend  on  the  line 

«  (  r  w 

termination  and  excitation.  The  terms  e  ’"•*  are  diagonal 
matrices  having  as  elements  e  -1,  where 

A  similar  development  for  the  line  voltage  (Vn(z.s))  can  be 
carried  out  to  determine  voltage  modes  and  a  propagation  equation  similar  to 
Equation  (6).  By  defining  a  characteristic  Impedance  matrix  as 


(Zc  ) 

n,m 


V) 


the  line  voltage  N-vector  can  be  txpresseo  using  the  same  constants 

(a*)  and  (cL“ )  as  In  Equation  (6): 
n  n 

(VnU.s))  -  (Zc  ){T 
n,a 

The  unknown  constants  (a*)  and  (5  )  are  determined  by 

n  n 

taking  Into  account  the  loads  at  each  end  of  the  line,  as  well  as 
the  excitation.  Consider  the  line  shown  In  Figure  2,  which  has 
lumped  voltage  and  current  sources  at  z  •  z$  ,  as  well  as  load  impe¬ 
dances  Ui  )  and  (Z-  )  at  z  ■  0  and  z  ■  I  respectively. 

I  p  p| 

On  the  section  of  the  lini  0  i  z  -  zs  Equations  (6)  and  (8) 
are  valid,  since  this  sect  _>n  of  the  line  is  source  free.  Similarly, 
for  i  z  i  l  similar  guatiois  are  valid,  but  with  different 
constants,  (dn)  .  By  relating  ,Vnu,s))  to  ( Xn(z ,s ) )  at 
z  ■  0  ,  and  z  •  t  through  the  load  impedance  matrices  and  by 
relating  the  discontinuities  of  (Vn(z,s))  and  (In(z,s))  to  the 
voltage  and  current  sources  at  z  ~  z$  ,  a  set  of  linear  equations 
can  be  developed  with  the  (in)  constants  for  section  of  line  as 
unknowns . 

Of  special  Interest  are  the  load  currents,  i.e  ,  (I  (OiS)) 
and  (in(t,s))  .  Using  the  solutions  for  the  (dR)  as  well  as 
Equation  (6)  for  z  ■  0  and  z  *  i  ,  the  load  currents  may  be 
expressed  as: 


/  -(>_  Jt  4  (Yn  JZ  .  \ 

...»  («  "•*  <<>  - «  •  <“„>)  <8> 
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where  the  terms  ( /  „(2S •  *))  and  (f^(zs*s))  represent  the  source 

terms  for  the  positive  and  negative  traveling  waves  on  the  multi- 

conductor  line.  Those  are  referred  to  as  combined  current  sources. 

since  they  have  the  dimension  of  current  but  arise  from  both  the 

applied  voltage  and  current  sources  at  2  »  z$  .  In  this  equation. 

the  terms  { f i  )  and  (f?  )  are  generalized  current  reflection 

'n.m  ‘•n.m 

coefficient  matrices  given  by 


<ri  > 

n.m 


k  1  *  izc  >r'  k  >  -  (2  )] 

s  n.m  n.m  *  l  n.m  n.m  J 


(10) 


for  the  load  at  z  *  0  ,  and  similarly  for  (T^  at  z  *  t 
with  (2?  )  as  the  load  Impedance.  As  defined  previously, 

(Iq  n)  Is  the  characteristic  Impedance  matrix  of  the  line. 

Notice  that  the  matrix  equation  in  Equation  (9)  has.  as  Its 
elements,  matrices.  Thus,  It  Is  referred  to  as  a  super  matrix  equation. 
The  double  dot  operator  (:)  Is  used  to  signify  the  product  between, 
two  super  matrices  by  first  treating  the  super  matrices  as  If  they 
were  regular  matrices  and  then  performing  matrix  multiplications  for 
each  of  the  Individual  multiplications  of  the  super  matrix  product. 

The  form  of  the  source  terms  In  Equation  (9)  can  be  shown  to  be 


<  * ;<*>> 


and 


<'n<*» 


(T  J-e 


1 

?  '  n.m 

♦  <! 


^n,m^2s 


(V^s)(zs.s)) 


(11) 


?  (Tr 

-  (I 


_)(*-z.) 

).e  s  .(T 


n.m 


((Zc 

V  n,« 


r1.(»i»)(*s.s)) 


(12) 


With  these  source  terms,  the  terminal  response  of  the  transmission 
line  can  be  determined  for  lumped  voltage  and  current  sources  at 
z  •  z$  .  For  field  excitation  of  the  transmission  line.  It  Is 
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necessary  to  consider  distributed  excitation,  as  opposed  to  the 

discrete  excitation  discussed  above.  This  can  be  regarded  as  a 

simple  extension  of  Equations  (11)  and  (12)  by  Integrating  over 

the  source  terns  (V^5h  and  (ij^).  Doing  this,  the  combined 
n  n 

current  sources  become 


(/>» 

fr>» 


i  /  (‘V.>'<V*K  j-'(»;(s)u.o) 


n,m 


♦  (i;(s)U.s))JdC 

-  (i;(s)  u,*)A  dc 


(13) 


(14) 


which  follows  directly  from  superposition.  Notice  that  now  the  voltage 
and  current  sources  are  per-unit*length  quantities,  and  hence  denoted 
by  a  prime.  These  quantities  must  be  determined  given  a  knowledge  of 
the  Incident  electromagnetic  fie  id  on  the  line,  as  well  as  a  knowledge 
of  the  transmission  line  cro^?- sectional  geometry.  This  Is  discussed 
In  the  next  section  of  this  report. 


IS 


SECTION  m 

DETERMINATION  OF  DISTRIBUTED  VOLTAGE  AMD  CURRENT  SOURCtS 


As  Indicated  In  the  previous  section,  the  terminal  (or  load) 
currents  of  a  multi conductor  transmission  line  can  be  evaluated  using 
Equations  (9),  (13)  and  (14)  If  the  distributed  voltage  and  current 
sources  (9^s)(z, s))  and  (I^s^(z,s))  are  known  everywhere  along 
the  line.  In  soaie  Instances,  such  as  a  small  aperture  or  other 
localized  source  close  to  the  transmission  line.  It  Is  possible  to 
approximate  the  solution  using  a  discrete  source  position,  as  In 
ref.  (9).  For  an  arbitrarily  Incident  plane  wave,  however,  this  Is 
not  possible.  Sources  distributed  over  the  entire  line  are  necessary. 

Consider  the  case  of  a  single  multiconductor  cable  In  free  space 
and  with  Impedance  terminations  at  each  end,  as  shown  In  Figure  3. 

Assume  that  In  this  bundle  there  are  n+1  wires,  with  the  n*ls* 
wire  being  the  reference  conductor.  The  electric  and  magnetic  fields 
In  the  vicinity  of  the  line  can  be  divided  Into  two  parts.  These 
are  the  Incident  components,  and  ff*nc  and  the  scattered 

components  C*  and  fi5  ,  such  that 

r-F'"0.!5  ns,) 

K  •  H*nc  *  S5  (15b) 

The  scattered  field  components  are  caused  entirely  by  the  Induced 
currents  and  charges  on  the  n*l  wires,  as  well  as  by  the  currents  on 
the  terminations.  The  scattered  fields  from  the  line  can  be  further 
subdivided  Into  three  different  classes.  There  are  TEN,  TE  and  TN 
transmission  line  modes,  which  are  produced  by  "transmission  line" 
currents,  having  the  property  that  the  components  of  the  total  current 
on  each  of  the  conductors  sum  to  zero. 

In  addition  to  these  currents,  there  are  "antenna  mode"  currents. 
These  are  currents  which  flew  on  each  wire  (but  with  a  different  magnitude 
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Figure  3.  Isolated  multi conductor  line  excited  by 
Incident  plane  wave. 


17 


for  each  wire.  In  general)  and  are  subject  to  the  constraint  that  the 
voltage  difference  between  any  two  conductors  In  a  transverse  plane 
Is  zero.  Furthermore,  these  currents  go  to  zero  at  the  ends  of  the 
line. 

Finally,  there  can  exist  guasl-statlc  current  and  charge  distri¬ 
butions  which  contribute  to  the  scattered  field  but  have  a  net  current 
or  charge  of  zero  on  each  conductor.  Although  these  latter  currents 
and  charges  do  not  play  a  role  In  computing  the  transmission- line 
response  directly,  they  are  Important  In  determining  the  coupling  of 
electromagnetic  fields  to  the  transmission  line. 

A  complete  and  rigorous  solution  for  the  field  Induced  currents 
on  the  multi  conductor  line  In  Figure  3  can  be  obtained  by  formulating 
and  solving  a  set  of  coupled  Integral  equations  for  the  wire  and  load 
currents,  given  a  particular  Incident  field.  In  many  cases,  however, 
such  a  complete  solution  for  t!-  current  Is  not  needed.  Fcr  lines 
which  are  long  compared  w^'.l:  the  wire  separation,  the  currents  due  to 
TE  and  TH  fields  attenuate  rapidly  from  the  loads  or  other  line 
terminations,  giving  rise.  theref."e,  to  a  current  distribution  which 
corresponds  primarily  to  the  TEH  -ents  plus  the  other  scattering 
currents  mentioned  above.  Moreover,  In  many  cases,  only  the  transmission- 
line  current  response  Is  desired  since  the  antenna  mode  currents  do 
not  contribute  to  the  load  response  In  the  general  case,  and  If  the 
transmission  line  Is  next  to  a  reference  ground  plane,  the  antenna 
mode  currents  are  not  excited  at  all.  Under  the  assim^tlon  that  the 
TE  and  TH  currents  are  negligible  and  neglecting  the  effects  of  load 
currents,  the  total  E  and  H  fKlds  In  the  vicinity  of  the  transmission 
line  can  be  written  as 

f  .  ,  fdnt  *  fEM  +  £st  (16a) 

and 

H  ■  H1nc  ♦  R*nt  ♦  R*0*  ♦  R*1  (16b) 
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where  the  subscript  (Inc)  refers  to  the  Incident  (or  free  space) 
fields,  (ant)  denotes  the  fields  produced  by  the  antenna  mode  currents, 
(TEH)  stands  for  the  fields  due  to  the  transmission- line  currents, 
and  (st)  Is  for  the  portion  of  the  fields  caused  by  the  static 
distribution  of  current  and  charge  on  the  wires,  determined  with  the 
condition  that  the  total  current  and  charge  be  zero  on  each  wire. 

Following  the  approach  used  in  ref.  (11)  for  single-wire  lines 
and  In  ref.  (7)  for  multi  conductor  lines,  Maxwell's  equations  can 
be  used  to  derive  a  v-1  relation  for  the  transmission  line  currents. 
Consider  a  uniform  section  of  multi conductor  line  shown  In  Figure  4. 

For  a  time  dependence  of  est  ,  Maxwell's  equation  may  be  written  as 

Vxt  -  -S  B  (17) 

and  on  a  path  .  from  the  reference  conductor  to  wire  1  (where 
dtj  represents  an  element  of  the  path,  and  fij  Is  the  normal  to  the 
path),  we  can  integrate  Equation  (17)  to  yield  the  following: 


This  result  Is  standard,  and  -;ts  derivation  will  not  be  repeated  here. 

Noting  that  the  line  Integral  of  the  electric  field  in  Equation 
(18)  Is  the  negative  of  the  voltage  between  the  two  conductors,  this 
equation  may  be  written  as 


11.  Lee,  K.S.H.,  "Two  Parallel  Terminated  Conductors  In  External  Fields, 
IEEE  Trans.  EMC.  Vol.  EMC-20,  No.  2,  pp.  288-296,  May  1978.  (A 
revision  of  ref.  6.) 


reference 

conductor 


n  *  dt  x  z 


Figure  4.  Cross  section  of  nu It 1 conductor  line  showing 
Integration  path  Cj  froai  point  a  to  b  . 
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As  discussed  by  Paul  In  ref.  (7),  the  term  Involving  the  B  field, 
which  arises  from  all  TEN  currents  on  the  multi  conductor  line  and  Is 
a  magnetic  flux  per  unit  length,  can  be  computed  In  terms  of  the 
Inductance  coefficient  matrix  elements  as 


( 


BTEM‘fi  dt  2  -  Ln  ^  ♦  L12  I2  ♦  ...  L,„  L  (20) 


In  ‘n 


where  Ij...ln  represent  the  currents  on  the  n  (non-reference) 
conductors. 

From  our  definition  of  the  "antenna  mode”  currents,  the  voltage 
between  wire  1  and  the  reference  Is  zero  for  these  currents,  which 
Implies  that  the  antenna  current  flux  term  Is  also  zero.  See  ref.  (12). 
Thus,  we  have  the  relation 


( 


B*nt*fl  dt  =  0 


(21) 


With  these  substitutions.  Equation  (19)  can  be  written  as 


-fiT'  -ju  a;,!,  » i;?i2 .  ...  i;nt„)  •  if  (Slnc  *  lst)-«  dt  (22) 


This  procedure  may  be  repeated  for  each  of  the  n  wires  In  the  bundle, 
and  the  resulting  equations  expressed  In  matrix  form  are 


12.  Frankel,  S.,  "Evaluation  of  Certain  Transmission-Line  Forcing 

Functions,"  AFWL-TR-78-T71 ,  Air  Force  Weapons  laboratory,  Klrtland 
Air  Force  Base,  NH,  (to  be  published). 


(23) 


d(vn) 

IT 


-s 


♦  s 


(i*nc 


♦ 


The  last  ter*  In  this  equation  has  dimensions  of  (volts/unit  length) 
and  Is  essentially  a  distributed  voltage  source  for  the  transmission 
line.  Denoting  this  by  (V^5^)  ,  we  then  have 


(24) 


where  the  relation  B  ■  u0H  has  been  used.  The  differentia'  equation 
for  voltage  and  current  In  Equation  (23)  then  becomes 


d<V 

dz 


♦ 


(v;(s)) 


(25) 


A  similar  manipulation  can  be  performed  using  the  other 
Maxwell  equation 


VxH  *  sc  £ 


(26) 


to  obtain  the  second  telegrapher’s  equation  containing  sources.  Applying 
this  to  the  contour  In  exactly  the  same  manner  as  In  ref.  (11), 
the  following  relation  may  be  derived. 
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CP 


By  Inserting  Equation  (16a)  Into  this  last  equation  and  noting  that 
the  antenna  mode  contributions  vanish,  since  by  definition  of  the 
antenna  currents,  ^r®n^*dt  *  0  and  f  H*n**fi  dt  *  0  ,  this 
equation  can  be  written  as 


♦  HTEM)-ft  dt 


■jif  <«,nc  * B 


or,  as  done  by  Lee  (ref.  11),  expressed  as 

-  £  f  C  io'-scC 

-4  *4 

Using  Equation  (20)  and  recognizing  that 
-Vj  ,  Eqautlon  (29)  becomes  "a 


(EinC+f*)-dt  (29) 
b  ^ 

E*dt  Is  the  voltage 


-  —  (l1  i 
Vi  dz  uiri 


L’lnln> 


»  -s 


(30) 


for  the  first  wire.  This  process  can  be  repeated  for  each  wire,  and 
the  following  matrix  equation  can  be  developed  for  the  transmission 
line  currents  (IR)  and  voltages  (Vn)  : 


1  d(in} 

i(L;.m)-dr-- -s  C<V 


Rearranging  terms  slightly  yields  the  second  telegrapher's 
equation 
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(32) 


d(M  (O 

— n~  ♦  S(C-  _).(v*)  *  (i;lSJ) 

<32  n  n  n 

where  the  source  term  (T^)  is  given  by 

<';,(5)>  ■  <r'"c  * 


(33) 


Note  that  in  deriving  this  relation,  the  assumption  that 


a*  j-(c: 


n,m 


n,m 


(34) 


has  been  employed,  a  result  which  implies  that  the  lines  are  within 
a  uniform,  homogeneous  dielectric  medium. 

In  an  Inhomogeneous  dielectric  region,  say  for  the  case  of 
each  conductor  having  a  separate  dielectric  jacket,  It  Is  known  that 
true  TEN  modes  cannot  exist.  However,  an  approximate  analysis  can  be 
carried  out  by  assuming  that  Equat  ins  (25)  and  (32)  are  applicable. 

The  validity  of  this  "quasl-TEM"  assumption  lies  In  the  reasonable 
comparison  of  theoretical  and  experimental  results  for  the  multi¬ 
conductor  system  (ref.  13). 

It  Is  to  be  noted  that  the  basic  telegrapher's  equations  derived 
here  for  the  transmission  line  currents  and  voltages  are  different  In 

form  than  those  developed  by  Paul  (ref.  7).  This  is  due  to  the  fact 

that  Paul  has  Integrated  from  the  center  of  one  conductor  to  the  other 

center,  not  from  one  surface  to  another  of  the  thin,  widely  spaced 

conductors  which  he  considers.  For  the  more  general  case  of  fat, 
closely  spaced  wires,  the  total  static  electric  and  magnetic  field  in 

13.  Chang,  S.K.,  F.M.  Tesche,  O.V.  Girl  and  T.K.  Liu,  "Transient 

Analysis  of  Multiconductor  Transmission-Line  Networks:  A  Comparison 
of  Experimental  and  Numerical  Results,”  AFWL-TR-78-152,  Air  Force 
Weapons  Laboratory,  Klrtland  AFB,  NM,  February  1979. 
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any  transverse  plane  must  be  used  to  compute  the  equivalent  line 
sources. 

Aside  from  a  difference  In  the  definition  of  the  unit  normal 
vector  h  ,  the  major  difference  between  the  formulation  of  lee  In 
ref.  (11)  and  the  present  analysis  Is  the  existence  of  an  additional 
antenna  mode  source  term  In  lee's  two-wire  analysis.  This  two-wire 
analysis  could  be  extended  to  a  multiwire  case,  and  thus  would  Imply 
the  existence  of  similar  source  terms  in  the  present  multi  conductor 
analysis.  As  discussed  by  Frankel  (ref.  12),  the  apparent  discrepancy 
arises  out  of  different  choices  for  the  "antenna  cu.-ent"  by  lee,  which 
thus  has  an  effect  on  the  remaining  transmission- line  current. 

As  stated  earlier,  our  choice  of  the  "antenna  current"  Is 
that  current  flowing  In  each  wire  which  produces  a  voltage  difference 
of  zero  between  any  conductor  and  another  at  any  transverse  plane  In 
the  line.  This  choice  is  also  used  by  Uchida  (ref.  14),  and  thus  leads 
to  a  decoupling  of  the  transmission  line  currents  from  the  antenna 
mode  currents. 

Although  explicit  expressions  for  the  voltage  and  current  sources 
have  been  developed  In  Equations  (24)  and  (33),  it  still  remains 
necessary  to  evaluate  the  scattered  static  fields  £*  and  H5  , 
before  the  source  terms  can  be  used  In  Equations  (13)  and  (14)  to 
determine  the  load  response  of  the  mul ti conductor  line.  To  determine 
these  source  terms.  It  Is  necessary  to  solve  two  static  boundary  value 
problems.  To  determine  the  current  source  In  Equation  (33),  it  Is 
necessary  to  solve  the  two-dimensional  static  problem  Illustrated  In 
rigure  5.  An  Incident  (free  space)  electric  field  strikes  a  collection 
of  conductors,  on  which  the  net  charges  are  zero.  A  static  scattered 
field  Is  produced  by  the  local  charges  Induced  on  each  wire,  and  the 


14.  Uchida,  H. ,  Fundamentals  of  Coupled  lines  and  Multiwire  Antennas. 
Sasaki  Publishing,  ltd.,  Sendai,  Japan,  1 967. 
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Reference  conductor 


Figure  5 


Cross  section  o*  multi conductor  cable  In 
Incident  F  field  showing  typical  field 
distribution  and  Integration  path  from  a 
to  b  .  Each  conductor  has  zero  net  charge 


integrals  in  Equation  (33)  are  then  evaluated  along  any  contour  fro* 
point  a  to  b  ,  using  the  total  scattered  field,  r*nc  ♦  E*  . 

The  solution  to  this  problem  for  the  multi conductor  case  is 
similar  to  the  two-wire  problem  discussed  by  Lee,  but  extended  to 
more  wires.  It  Is  solved  by  looking  for  the  solution  to 

V2*  «  0  (34) 


exterior  to  the  wires,  with  the  condition  that  ♦  «  constant  on 
each  of  the  conductors  subject  to  the  constraint  that 


on  each  conductor,  1  ,  and  that  at  infinity,  the  potential  Is 

*  -r*E  (36) 

H*re  $^nc  represents  the  Incident  or  free-space  potential  field  in 
the  absence  of  the  transmission  line.  Once  this  equation  is  solved 
(usually  by  numerical  means)  the  potentials  of  each  wire,  ,  can  be 
determined,  and  the  integrals  of  Equation  (33)  can  be  determined  directly 


(E 


Inc 


E*)' 


dt. 


-u, 


Vi* 


(37) 


It  is  possible,  however,  to  express  the  integral  In  Equation  (37) 
in  a  simpler  form,  using  only  the  Incident  field,  ,  and  a  vector 

equivalent  distance,  h^  ,  in  a  manner  similar  to  that  of  ref.  (11). 
Consider  an  auxllliary  problem  which  has  a  potential  field  given  by 
and  is  defined  by  the  relations 


2  * 

vU 


0 


(38) 


wi th  a* 


constant  (but  unknown)  on  each  of  the 


conductors  of 
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the  multi conductor  bundle,  and  with 


/  V*' 


for  all  conductors  except  for  the  1tn  conductor  and  the  reference 
conductor,  where  we  have  the  constraint 


J  "5T  dst  ‘  ‘  T 


wire  1 


/  ? 

reference 


U  •>+ 

ZClidS  ■  -i 
*Vl  "4l  f 


The  solution  to  this  auxiliary  problem  can  be  used  to  find  the 
field  excitation  of  the  transmission  line  by  using  Green's  Identity, 

♦•V2*  -  ♦V2**  «  0  (42) 

and  applying  Gauss'  theorem  to  give  the  expression 


all  conductors 


/  U  -  ♦  -  *r ) dS  *  /  (**  U  ■  ♦  ^5f] 


dS  •  0 


where  Is  a  closed  surface  at  Infinity.  Using  the  facts  that 
♦  and  e*  are  constant  on  the  conductors,  that 


(Hi 

j  >"j 


dS  -  0 


for  all  conductors  except  the  1tn  and  the  reference  conductor,  and 


/(♦*&- ♦£)'1S' •  /  (•"K  n£)dS  •  - 1  J  ♦"*<’*“ 


wires 


wires 
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where 


ia 

3n 


(46) 


Is  the  charge  density  on  each  conductor  for  the  auxiliary  problem  and  Is 
a  known  quantity.  Equation  (43)  can  then  be  expressed  as 


Using  Equation  (40)  and  the  relation  #1nc  *  -E*r  ,  this  last 
equation  takes  the  form 


<‘l  -  W 


where  the  vector  h^  is  defined  as 

i  roj  dSj  ♦  i  r  dS?  *  ...  ^ 
*S,  i,  S_ 


r  "JUl  dSn.l 


/ 


(48) 


(49) 


With  this  expression.  Equation  (37)  can  be  conveniently  expressed  as 


Vi 

ji  (f1nc  ♦  E^-dT,  •  rinc*hi  (50) 

a 

and  the  N  vector  equivalent  current  source  becomes 

<5,> 

The  vectors  h^  are  referred  to  as  the  "field  coupling  vectors" 
for  the  line,  and  also  as  the  "effective  height"  of  the  conductors. 
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Physically,  they  correspond  to  the  vector  distance  between  the  charge 
centroids  on  the  multi  conductor  system,  given  a  total  charge  Q  on  the 
1th  conductor,  -Q  on  the  reference  conductor  and  zero  net  charge 
on  all  others.  Figure  6  illustrates  these  relationships. 

For  the  case  of  thin,  widely  separated  wires,  the  vectors 
h^  are  simply  the  distances  from  the  center  of  the  reference  con¬ 
ductor  to  each  of  the  wires'  centers.  For  more  closely  spaced  wires, 
the  field  coupling  parameters  must  be  calculated,  using  the  Integral 
equation  approach  outlined  by  Girl  in  ref.  (15). 

A  similar  procedure  can  be  carried  out  for  determining  the 
distributed  voltage  source  in  Equation  (24)  by  solving  a  magnetostatic 
problem.  The  details  of  this  are  identical  to  that  described  by  Lee 
(ref.  11),  modified  by  the  presence  of  more  than  just  two  conductors. 
The  results  are  that  the  same  field  coupling  parameters,  h^  ,  that 
are  used  for  the  electric  field  calculations  may  be  used  for  the 
magnetic  fields.  This  results  In  the  following  equation  for  the 
distributed  voltage  source. 

(Vn(S))  1  Sko  (<VV‘*1nC)  («) 

The  preceding  discussion  has  been  for  the  field  excitation  of 
an  Isolated  multiconductor  line,  in  which  one  of  the  conductors  In 
the  bundle  serves  as  the  reference.  An  often  encountered  situation, 
however,  is  not  this  configuration,  but  one  with  an  n-wlre  bundle 
next  to  a  flat,  conducting  ground  plane.  For  this  case,  the  ground 
plane  serves  as  the  reference  conductor,  and  the  antenna  mode  currents 
aie  not  excited. 


15.  Girl,  D.V.,  F.H.  Tesche,  and  S.K.  Chang,  "A  Note  on  Transverse 
Distributions  of  Surface  Charge  Densities  on  Multi  conductor 
Transmission  Lines,"  AFWL  Interaction  Notes,  Note  337,  April  1,  1978. 
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Figure  6.  Cross  section  of  Isolated  n+1  wire 

multiconductor  line,  showing  field  coupling 
vector  for  the  1th  conductor. 


For  this  case,  the  field  coupling  parameters  are  still 
calculated  as  above.  For  example,  as  show  In  Figure  7,  the 
coupling  parameter  h^  Is  calculated  by  placing  a  charge  Q  on 
Mire  1  and  no  net  charge  on  the  other  wires.  By  Image  theory, 
there  Is  an  image  charge  of  -Q  on  the  image  of  wire  1  and  the 
resulting  charge  centroids  may  be  computed.  The  coupling  parameter 
vector  Is  directed  away  from  the  ground  plane  and  has  a  magnitude  equal 
to  the  shortest  distance  from  the  ground  plane  to  the  1th  wire's 
charge  center. 

In  this  case,  note  that  the  incident  fields  f*nc  and  H*nc 
which  are  used  In  Equations  (51)  and  (52)  must  include  the  reflection 
effects  of  the  ground  plane.  Thus,  If  f^nc  and  Hinc  represent 
the  free  space  fields  In  the  absence  of  the  ground  plane,  the 
exciting  fields  of  the  line  to  be  used  In  the  above  equations  are 


En  -  2(E,nc-n) 
and 

Ht  «  2(t*n)-Hinc 


(53) 

(54) 


where  fi  Is  a  unit  normal  to  the  plane,  6  Is  the  direction  of 
propagation  of  the  Incident  wave  and  the  subscripts  n  and  t 
represent  field  components  normal  to  and  parallel  to  the  ground 
plane,  respectively. 
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SECTION  IV 

EXCITATION  FIELDS  DUE  TO  INCIDENT  PLANE  HAVE 

The  expressions  for  the  distributed  current  and  voltage  sources 
In  Equations  (51)  and  (52)  In  the  previous  section  are  quite  general 
and  depend  only  on  the  local  Incident  electric  and  magnetic  fields 
on  the  line.  One  type  of  Incident  field  which  is  useful  to  consider 
Is  a  plane  wave  of  arbitrary  angle  of  Incidence. 

Consider  a  single  transmission- line  tube  being  Illuminated  by 
a  plane  electromagnetic  field.  As  shown  In  Figure  8a,  the  tube  Is  In 
the  i  direction  and  the  k  vector  of  the  Incident  field  arrives  with 
angles  with  respect  to  the  r  axis  and  PQ  ,  which  Is  the  Incli¬ 
nation  angle  of  the  Incident  field.  Two  different  polarizations  of 
the  Incident  field  are  possible,  and  are  denoted  as  TE  and  TM, 
respectively.  The  TE  case  occurs  when  the  Incident  E  field  Is 
perpendicular  to  the  plane  of  Incidence,  which  Is  defined  as  the 
plane  formed  by  the  k  vector  and  Its  two-dimensional  projection  In 
the  x-y  plane.  The  TH  case,  conversely,  occurs  when  the  E 
field  lies  within  the  plane  of  incidence.  Figure  8b  Illustrates 
these  different  polarizations. 

For  both  of  these  polarizations,  the  field  components  at  the 
multi  conductor  tube  can  be  expressed  as  follows: 

TE  Fields 

H1nc  ,  _H1nc  s1n  ^  E1nc  .  _E1nc  $if)  ^  CQS  ^ 

H1nc  .  H1nc  CM  ^  £1nc  .  _E1nc  s1(|  ^  $in  ^ 

H*nc  -  0  E*nc  •  Einc  cos  0O 
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Incident  propagation  direction 


y 


(TM  Polarization) 


(b) 


Figure  8.  Geometry  and  polarization  of  the 
Incident  plane  wave. 
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As  seen  from  Equations  (51)  and  (52),  the  important  quantities 
for  determining  the  distributed  sources  are  the  electric  field 
component  parallel  to  the  vectors  h^  ,  and  the  magnetic  field 
component  perpendicular  to  h,  .  Consider  the  geometry  shown  in 
Figure  9.  The  1  conductor  Is  shown  with  its  coupling  vector  having 
an  angle  9^  with  respect  to  the  chosen  x  axis,  and  a  magnitude 
h^  .  For  this  case,  the  components  of  the  electric  field  in  the 
direction  parallel  to  h^  are  given  by  the  following  expressions 
for  the  1th  conductor: 


and 


E*nC  sin  rtj  ♦  Ex  cos  0^ 

4  nr 

E  (sin  9.  cos  6  -  cos  ft.  sin  ft  sin  ) 

<0  too 

(TE  polarization) 

-E1nc  cos  ftj  cos  (TW  polarization) 


*H  cos  ♦  Hx  sin 

H1nc  sin  9^  cos  (TE  polarization) 

-H^nc  (cos  9,  cos  9„  *  sin  ft.  sin  0^  sin 

io  too 

(TM  polarization) 


(53a) 

(53b) 


(54a) 


(54b) 


With  these  field  components,  the  distributed  vector  current 
and  voltage  sources  take  the  form 
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reference 


Figure  9.  Cros*  -i.  on  of  multi  conductor  line  showing 
field  coupling  parameter  and  pertinent 
field  components. 


I 


and 


(i;(s,)--s  (c;  )  (k  t„  ) 

'  n 

•  s%  <\ ",  ) 
n 


(55) 

(56) 


and  should  be  used  In  Equations  (13)  and  (14)  to  evaluate  the  line 
response  for  an  Incident  plane  wave. 
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SECTION  V 
CONCLUSIONS 

This  report  has  presented  a  discussion  of  the  field  excitation 
of  multi  conductor  transmission  lines.  First,  a  general  expression 
for  the  current  response  at  the  terminations  ofanN-wIre  multicon¬ 
ductor  cable  has  been  developed  In  terms  of  distributed  voltage 
and  current  sources.  In  Section  III  relationships  between  these 
sources  and  the  total  static  electric  and  magnetic  fields  In  the 
vicinity  of  the  transmission  line  are  then  derived.  These  are  then 
related  to  the  free  space  (or  Incident)  fields  through  a  vector  field 
coupling  parameter  or  equivalent  separation  of  the  lines.  Finally, 
Section  IV  expresses  the  distributed  source  terms  for  the  multi¬ 
conductor  line  In  terms  of  the  angles  of  Incidence  and  polarization 
of  an  Incident  plane  wave. 

This  work  expands  upon  the  past  studies  of  field  excitation 
of  two-wire  transmission  lines.  The  field  coupling  parameters  for  a 
multiconductor  line  are  seen  to  be  determined  from  a  series  of  cal¬ 
culations  Involving  specifying  a  zero  net  charge  on  all  conductors 
except  the  reference  and  the  conductor  for  which  the  coupling  para¬ 
meter  Is  being  determined.  It  Is  noted,  furthermore,  that  the 
excitation  of  the  line  depends  strongly  on  the  line's  orientation 
with  respect  to  the  Incident  fields. 
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